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ABSTRACT 

Evidence of small-scale condensations in the magnetic molecular clouds has been ac¬ 
cumulating over the past decades through radio and optical/ultraviolet observations. 
The origin and shape of these small-scale condensations is a disputable issue. Nejad- 
Asghar & Ghanbari (2004 hereafter NG) have recently studied the effect of the linear 
thermal instability on the formation of fluctuations in molecular clouds. The authors 
inferred that under certain conditions (e.g., depending on expansion or contraction of 
the background) thermal instability and ambipolar diffusion can produce spherical, 
oblate, or prolate condensations. 
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1 INTRODUCTION 


Firstly, about presence of these fluctuations. We have two 
method to find them: (1) direct imaging of nearby clouds re¬ 
veals substructures on scales to lengths of 0.0Ipc and masses 
of O.OIMq (Peng et al.l998; Sakamoto & Sunada2003), (2) 
studies of the time variability of absorption lines indicates 
the presence of fluctuations on scales of 10“^pc (5-50 AU) 
and masses of IO^^Mq (Moore & Marscherl995; Boisse et 
al.2005). 


How these fluctuations may be formed and is the ther¬ 
mal instability important? To answer the question, first, we 
must consider time-scales. According to the prior results of 
Gilden (1984), thermal instability time-scale in molecular 
clouds is in the order of 10® — 10"^ year. On the other hand, 
Larson (1981) showed that dynamical time-scale for tran¬ 
sient structure of molecular clouds is in order of 10^ year. 
Thermal instability time-scale is very smaller than dynam¬ 
ical time-scale. Thus, it is reasonable to consider thermal 
instability as an important mechanism in formation of small 
fluctuations in molecular clouds. 


It is now accepted that the magnetic field is very im¬ 
portant in dynamical evolution of all ISM. It will affect on 
ions directly and on neutral particles indirectly, via colli¬ 
sion with ions. In molecular clouds, ionization degree is very 
small, thus, drift of neutral particles between tied ions is im¬ 
portant. It can heat the medium. In this way, for complete 
investigation of the thermal instability in molecular clouds, 
we must consider plasma drift (ambipolar diffusion). 


2 NET COOLING FUNCTION AND 
PERTURBATION 


The most important parameter in thermal instability is 
the net cooling function. A good fitted result for cooling 
in molecular clouds is A(p, T) — Aop^T^ (Goldsmith & 
Langerl978; Neufeld et al.l995), where the values of S and 
ff is given by the Fig. 1 of NG. 

About heating rate, we turn our attention to gravita¬ 
tional expantion/contraction work and ambipolar diffusion 
heating rate. The gravitational heating rate is found by set¬ 
ting the rate of compressional/expansional work per parti¬ 
cle, equal to the rate of change of gravitational energy per 
particle. For ambipolar diffusion heating rate, we have 


Fad = (1) 

pn 

where = rjep^'^'' is the drag force (per unit volume) with 
rj as collision drag, and Vd is the drift velocity of ions relative 
to neural particles. In order of magnitude, if kBq changes on 
a typical scale of A, and if we choose pi ~ ep)) (Umebayashi 
& Nakanol980), ambipolar diffusion heating rate, in a good 
estimate, is given by 


iAD =IoP ? 


IGrr^peA^ 


( 2 ) 


In this way, the net cooling function in molecular clouds is 
given by (for more detail see NG) 


H(p,r) = AopV'® - Fop®/^ - rop“(^+"^ (3) 

Instability condition depends on the different values of /?, 5, 
and V. 

We use the perturbation method to investigate the oc- 
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Figure 1. Regions of stability(x), spherical instability(•), and 
oblate instability(—) for the case of Gg = 0.1 and Eg, = 0, for a 
typical value of o = 5.0 and D = 1.0. 


Figure 2. Regions of stability( x), spherical instability(•), and 
oblate instability(—) in the expanding background, for a typical 
value of o = 5.0 and D = 1.0, with (a) Eg = 0.01 and Gg = 0.1, 
and (b) Eg = 0.1 and Gg = 0. 


currence of thermal instability in molecular clouds and for¬ 
mation of fluctuations. We choose an expanding/contracting 
molecular cloud with uniform magnetic field Bq. Perturba¬ 
tion on this medium results to a five-degree linear character¬ 
istic equation, which its coefficients depend on the impor¬ 
tance of ambipolar diffusion and gravitational heating rate 




Pad _ P grav 

A(no,2o)’ ^ A(no,To)’ 


(4) 


and some time-scales that are investigated in details at sub¬ 
section 3.2 of NG. By introducing the non-dimensional quan¬ 
tities 


V = hTs, CTp 
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(5) 



a = {A^f,D = 

TAL 


Gg = (^)^ Eg = li, (6) 

Tg Tg 


we use the Laguerre method to find five roots of the charac¬ 
teristic equation (?). 


Figure 3. Regions of spherical instability(•), and prolate 
instability(l) in the contracting background, for a typical value 
of o = 5.0 and D = 1.0, with (a) Eg = —0.01 and Gg = 0.1, and 
(b) Eg = —0.1 and Gg = 0. 


3 RESULTS AND PROSPECTS 

If we consider the effect of the self gravity without expan¬ 
sion/contraction of the background, the isobaric instability 
criterion (line OA of Fig. 1) is modified, because, self-gravity 
causes to increase the internal pressure. 

If the background is expanding {Eg > 0), its expansion 
energy causes to stabilize the medium. This case is shown 
in Fig. 2 for a typical value of a = 5.0 and D = 1.0. In the 
isentropic instability criterion (line OB), expansion energy 
causes to stabilize the medium in the direction of the mag¬ 
netic field and perpendicular to it. In the isobaric instability 
criterion (line OA), it only causes to stabilize the medium 
in perpendicular to the magnetic field, corresponding to de¬ 
creased pressure via ion-neutral friction. 

For contracting background {Eg < 0), contraction en¬ 
ergy injected to the medium, thus, its stability is decreased 
and converted to a prolate instability. Diffusion of neutrals 
relative to the freezed ions in the perpendicular direction of 
the magnetic field is the reason of this prolate instability. 
This case is shown in Fig. 3 for a typical value of a = 5.0 
and D = 1.0. When the parameters of a magnetic molecular 


cloud set, locally, in this region of or — Up plane, prolate 
condensation may be produced via thermal instability. 

Investigation of this problem in the linear prospect is 
not completely OK, because, ambipolar diffusion is a non¬ 
linear effect and density fluctuation ratios are in the order 
of 10. But, we see that linear study can give us a good idea 
about the conditions for occurrence of thermal instability. 
For complete study we must work in nonlinear regime. 
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